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Abstract 

Weakly activated signaling cascades can be modeled as linear systems. The input-to- 
output transfer function and the internal gain of a linear system, provide natural measures 
for the propagation of the input signal down the cascade and for the characterization of 
the final outcome. The most efficient design of a cascade for generating sharp signals, is 
obtained by choosing all the off rates equal, and a "universal" finite optimal length. 

1 Introduction 

Protein kinase cascades are major functional modules used by cells to translate signals gener- 
ated by receptor activation into diverse biochemical and physiological responses ^H^. Highly 
conserved throughout evolution and across species, the kinase cascade motif participates in the 
control of many processes, including cell cycle regulation, gene expression, cellular metabolism, 
stress responses, and T cell activation. For this reason, control of kinase cascades by therapeu- 
tic intervention has become an attractive area for drug discovery, particularly in the areas of 
cancer and inflammation |§] . 

Some four mitogen-activated protein kinase (MAPK) signaling cascades have been found in 
yeast ^7j and at least a dozen MAPK cascades have been identified in mammalian cells |12| . 
The intensive study of MAPK pathways has prompted efforts to characterize these systems 
theoretically (see, inter alia, E3 HJ H3 H3 E])- in this paper, we will utilize 

concepts and methods from the theory of linear control systems to characterize kinase signaling 
cascades, and in particular the MAPK pathway, in order to understand how the number of 
kinases in a cascade and their individual enzymatic activities can affect the pathway in its role 
cLS cl SI gnal transducing module. 

Let R denote the input signal, Xj the inactive (nonphosphorylated) form of kinase i and 
Xi the active (phosphorylated) form of kinase i. The rate constant (or "on" rate) for the 
z-th kinase phosphorylation will be denoted by &i, and the dephosphorylation rate constants 
(or "off" rate) will be denoted The input signal R might represent, for example, the 
concentration of activated receptors, and the dynamics of the signal transduction pathway may 
be modeled as follows (see [TU]): 

^ = a 1 RX 1 -0 1 X 1 , ^=a i X i . 1 X i -/3 i X i , i = 2,...,n. (1) 
at dt 

Assuming that the total amount of kinase i remains constant, that is, X \ + Xi = X tot; j the 
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differential equations (pQ) can be rewritten as 



and 



' /A ' Oi^ifl-^-) - i,A,. ,-■> „. (2) 

L tot,i 



dt V X 



where aj = diX tot i- Throughout this paper we will focus on the case of weakly activated 
pathways, by which we mean a low level of kinase phosphorylation, that is 

Xi<X totti (3) 

In this case the equations (J2J) are simplified to a linear system of the form: 

^ = a 1 R-(3 1 X 1 , ^ = aiXi-! - PiXi, % = %... n. (4) 

In Section [21 we will describe how to compute the transfer function and internal gain for this 
system and then in Section |3] we will define a set of measures for the output signal, which closely 
follow those discussed in ^0]- in Section H] we prove that the most efficient cascade design, for 
generating sharp signals, has equal on rates and a finite length depending only on the cascade's 
internal gain. In Section positive feedback from the last activated kinase to the first is added 
to the cascade, and the optimal design is re-examined in this new context. Finally, in Sections^ 
and we briefly discuss the effect of delays along the cascade and how to check the cascade's 
stability to random small perturbations. 



2 The input-to-output transfer function 

We will consider the signaling cascade (jU as a system with an input R, and an output which 
will be some function of the concentration of the last kinase X n . Specifically, the output will 
be the "effective" integral of X n , or in other words, the cascade will be extended one more step 
to include a "leaky" integrator: 

dX n+i _ y fY 

where the output is X n+ \. The variable X n+ \ expresses the effective concentration of the last 
kinase (minus losses due to degradation or inactivation of X n , for instance). Note that the case 
£ = recovers X n+1 = /* X n (t') dt'. 

The model for a weakly activated signal transduction cascade may then be written in the 
more compact form, 

jy 

(t) = AX(t) + BR(t), Y(t) = CX(t), (5) 

dt 
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input= R 




X n+I =output 
Figure 1: A model of a MAPK cascade. 



where X = (X\ , X2 , . . . 

A e R(«+i)x(n+i) B G 



X n , X n+ i)' is a column vector whose elements are always nonzero, and 
j(n+l)xl and C e K ix(n+i) are the matrices 
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and 



C = (0 



1 



(7) 



It is well known (see or an y other book on control systems) that, for a system such 

as (J5J), the output can be computed directly as the convolution between the input signal R 
and the impulse response of the system. The impulse response of the system is the output 
corresponding to a single input pulse. If we let G denote the impulse response (and assuming 
that the system starts at rest, with initial condition X(0) = 0), then 

Y(t) = (G*R) (t). 

The impulse response, G, characterizes the action that the internal structure of the system will 
have on any input, such as the filtering of certain frequency components, and the amplification 
or dampening of the signal. Biological inputs may take many different forms, such as single 
pulses, slowly decaying signals, constant stimuli applied for a certain time interval, or oscillatory 
signals. Thus, it is appropriate to have a model in which the output signal is obtained as a 
convolution of the input R (which may take many forms) and the transfer function G (which 
depends only on the intrinsic kinase activity parameters, and needs to be computed only once). 
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A very convenient way to analyze system © is to convert it to the frequency domain, by 
application of the Laplace transform operator. The Laplace transform of the impulse response is 
called the transfer function of the system, and it provides a simple linear relationship between 
the Laplace transforms of the input and the output, as well as also providing a measure of 
amplification/dampening of the input signal. The transfer function is given by a simple formula 
in terms of the matrices A, B and C as summarized in Appendix [DJ For this cascade system, 
we will carry out the Laplace transforms in detail so as to gain some insight into the internal 
structure of the system. 

The Laplace transform of X will be denoted by X, and is defined as 



Ms) 



e~ st Xi(t)dt, and R{s) 



-st 



R(t)dt, 



where s is a complex number s = s re + jlu (j is the imaginary number v— T) and takes values 
in an appropriate region of convergence. 

Applying the Laplace transform operator to both sides of equations (JIJ), assuming that 
X(0) = 0, and recalling the properties of the Laplace transform (see Appendix[Dj), we have: 

sX^s) = aiR(s) -piX y {s) 
sXi(s) = aiXi_i(s) - PiXi(s), i = 2,...,n 
sX n+1 (s) = X n (s) - £X n+1 (s), 



which yields 



and 



Xi(s) 



R(s), Xi(s) 



Xj_i(s), i = 2, ...,n, 



X n+ i(s) 



s + 



X n (s). 



In this way, we may view the cascade as a sequence of n steps, the output of the step i — 1 
becoming the input to step i. 
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Figure 2: Transfer functions at each step. 



For each single step in the cascade, the input is X{-\ and the output is X{, and they are 
related by a multiplicative factor, which is in fact the transfer function for the step i: 



Gi(s) 



CX; 



s + A 



i = 2, ... ,n. 



For the whole cascade, the input is R and the output is X n+ \, and it is easy to see that the 
transfer function for the total system is the product of the transfer functions at each step: 



G{s)=G 1 {s)---G n+1 {s) 



1 



ol\ ■ ■ ■ a n 



s + £(s + ^)---(s + /3 n y 



(8) 
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Therefore 



Y(s) = X n+1 (s) 



R(s), 



and the actual output may now be obtained by the inverse Laplace transform. Alternatively, 
even without knowing the exact form of the output, that is, the function Y(t), it is still possible 
to further characterize the properties of the system, through the 2-norm of the functions Y and 
R. Define the 2-norm of the function Y and the 2-norm of the Laplace transform Y by 



1*1 



\Y(t)\ 2 dt 



and ll^lb 



1 

2^ 



\Y{ju)\ 2 duo 



and similar expressions hold for R and R. (Note: from now on we will assume that the signals 
are defined only for positive times, that is, Y(t) = for t < 0). The 2-norm ||i^||2 provides 
a measure of the strength of the signals (in analogy to the energy of a mechanical system). 
Indeed, these norms provide a very convenient way to relate the input and output because, 
from Parseval's Theorem, the 2-norm of a function, equals the 2-norm of its Laplace transform, 
and therefore 



\Y\\ 



111 



l-RI 



\R\\ 



without the need to compute inverse transforms (a very helpful fact, since in general the inverse 
transforms may not be simple to compute). 

Another useful measure is the infinity norm of the transfer function, that selects the least 
upper bound of the absolute value of G, 



\G\ 



sup|G(ju;)|. 



As shown in Appendix [Dl a very useful estimate for characterizing the relative strength of the 
input and output signals is: 



\Yh < \\G\\ 



I Jill 



(9) 



where it is immediately apparent that the infinity norm of the transfer function gives an upper 
bound for the amplification of the input signal throughout the cascade. Moreover, the infinity 
norm ||G||oo is in fact the smallest number that satisfies Q, for all input /output pairs (that is, 
pairs (R,Y), where Y is the output corresponding to the input R). 

To compute the infinity norm of the transfer function for the whole cascade, note that 



\Gi(ju)\' 



ad 



a 



\jv + Pi 



" 2 + P?- ft 



< -4, for all u G (— oo, oo) 



and the equality holds for uj = 0. Therefore 



I CI 



1 ai • • • a n 
£lh---Pn' 



(10) 



A necessary condition for amplification of the signal to occur is that ||G||oo > 1. Moreover, since 
t is essentially an independent parameter, introduced for the purpose of defining a reasonable 
measure of the output, we can say that amplification of the input signal occurs only if 



ai ■ ■ ■ a n > fli ■ ■ ■ p n . 



(11) 
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Recall that a, = djX totj j, where X tot j is the total concentration of the ith kinase and on is the 
(true) rate of phosphorylation. Therefore, we still expect that oti < fy, i = 1, . . . , n, as should 
be the case for a weakly activated pathway. 

The norm HGHoo is often called the internal gain of the system which, through expression @, 
provides a useful and easy way to compute the input-to-output strength relation. For example, 
if a MAPK cascade has a "5-fold amplification", then its internal gain is HGHoo = 5. 

Note that, in the case where I = 0, the internal gain HGHoo is infinite — meaning that, in 
at least one step ( X n — » X n+ \ ) there is no degradation term. Then the estimate © contains 
no useful information. However, for t = 0, we have 



Y(t) = X n+1 (t) = / X n (t') dt', 
Jo 

and we also have an estimate for the "strength" of the signal X n , since 

Oil " " " a n 



\X. 



n || 2 
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R 
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3 Signaling time, signal duration and signal amplitude 

Some basic quantities which serve to characterize a signal transduction system are: the overall 
amplification from the input to the ouput; the duration of the output signal; and the time 
it takes the input signal to traverse the cascade. There are several possible definitions and 
estimates of these quantities: here we extend the definitions given by ^Uj, embedding them in 
the context of frequency-domain analysis, and generalizing them to arbitrary inputs. 

To be concise, let us identify the cascade (jHJ by its parameters, and associate with it the 
following (2n + l)-tuple: 

C := (n,«i, ... ,a n , . ,/3 n ), 

where it is assumed that n E N and oti, Pi are positive real numbers, for i = 1, . . . , n. We will 
also introduce the notation IA for denoting the set of inputs. 

Definition 3.1 For system (j^J), with parameters C and a leak factor I > 0, for each input R, 
the signaling time, r, and the output signal duration, a, are given by 



. ^ „ . d In Y , . 

r(C,i,R) := - —^~(s) 



o(C,£,R) :-- 



J s=0 



d 2 my 

ds 2 



J s=0 



The signaling time to step i and the signal duration at step i, i < n are given by: 



„n d In Xi , . 

n(C,R) :=__ ^_ i( a ) 



o-i{C,R) 



s=0 



d? hxXj 
ds 2 



J s=0 



To understand the significance of these definitions, recall the properties of the Laplace transform 
and compute (with Y(t) = for t < 0): 



Y(0) 



Y(t)dt, ^(0) 
ds 



d 2 Y 

tY(t)dt, ^(0) 



t 2 Y(t)dt 
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and thus we recover expressions (4) and (5) of reference ^U] 



_ f~tY(t)dt 2 



f °°Y(t)dt \fi°Y(t)dt) 



where r can be regarded as the expected value (of the time to traverse the pathway), and a as 
the corresponding variance. 

An estimate of the amplitude of the output signal, as given in equation (6) of reference [TUj . 
is the value S, such that S x 2a = J °° Y(t)dt. Again we propose a more generalized notion, 
suggested by the input-to-output estimate ©, that takes advantage of the easily computed 
internal gain of the system, and also incorporates the strength of the signal. 

Definition 3.2 For system (jSJ), with parameters C and a leak factor £ > 0, for each input R, 
the signal amplitude is given by 



A(C,£, R) 



IIGelloo 
a{C,£, R) 



\R\\ 



(12) 



where Gq is the transfer function (JHJ). 

A may also be regarded as the amplitude of a constant signal of duration a, but Definition ^. 21 
differs from the definition of amplitude given in ^U] in essentially three points: 

1. the meaningful quantity for measuring the amplitude is not the integral jY(t) dt (which 



computes the area under the curve Y(t)), but rather the 2-norm y f \Y(t)\ 2 dt, which 
computes the strength of the signal; 

the amplitude is proportional to the product of the internal gain of the system, and 
the 2-norm of the input. This simplifies calculations since, for each cascade, the \\G\\ is 
computed only once and ||i?||2 is computed for each input signal; 

the product ||G||oo||-R||2 is used as an estimate for ||y||2, but we know (see Appendix) that 
1 1 G | |oo is the least factor that satisfies the inequality ||y||2 < K ||-R||2- In fact, jS] shows 
how to construct examples of inputs for which the equality is approximated. For instance, 
for any e > 0, the input depicted in Figure |3 



R(t) = 2— sinet, 



with r 



for t > 0, 



(13) 



has unit norm, i.e., ||-R||2 
in the Appendix. 



1, and satisfies ||y||2 ~ HGHoo, for e small enough, as shown 



We remark that these definitions are valid not only for the special case when A, B and C 
are of the form specified in equations . (|7 |) . but in fact they are valid for any linear system of 
the form ©. For instance, in Section [51 we compute these quantities for the case when there is 
positive feedback from the last to the first kinase. We next explicitly compute these quantities 
for the special case when A, B and C are of the form (jEJ) and 0, and £ = 0: 



t(C,£,R) 



1 1 



d ]nR 

ds 



(14) 
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Figure 3: An input that satisfies ||-R||2 = 1 and ||V||2 ~ H^Hoo, with e = 0.2. 



a(C,£,R) 



1 - 1 

m p. + J2 -p? + q( r )> where i( R ) 



d 2 lnfl 

ds 2 



J s=0 



(15) 



A(C,£,R) 



\R\\ 



(16) 



In the case £ = 0, the quantities r, a and A may be computed for Y = X n . The expressions 
are very similar, except that all the terms in I vanish. 

Example 3.3 A typical input is a decaying exponential R(t) = Roe~ xt , with 
Rq a, . Rq d \nR 



R\ 



R(s) 



(0) 



1 



q(R) 



2A' ^" y "' s + y ds 
A "peak" -like input may be represented by R(t) = Rote~ xt , with 

R Q d InR 



1 

A2' 



\R\\ 



(0) = "f. ?(*) = JS- 



Cs + A) 2 ' ds 

For a constant signal, of magnitude Rq, applied for an interval of time Tq, we have: 



\R\\2 — Roy Tq, R(s 



Ro 



1 



-sT 



d InR 

ds 



(0) 
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4 Cascade design optimization 



From the analysis of the quantities r, a and A, defined in Section |3J we can explore the 
signaling efficiency of kinase cascades. The definition of an "efficient" response may depend on 
the particular biological context, but it typically involves the relationship between the length 
of the cascade, the amplitude of the signal and its duration. A question posed in ^0] is 
whether cascades can respond with sharp signals, i.e., simultaneously of short duration and 
high amplitude. Our model provides a definite answer to this question. 

As we have seen, our linear model has a gain that depends on the length of the cascade 
and the values of the on/off rate constants, but doesn't depend on the input. As a starting 
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point, we may think of the family of cascades that have the same value for the internal gain, 
say K, and examine their length, the distribution of the "on/off" rates and signal amplitude 
and duration. The problem we would like to study is then: 

(P) For each fixed internal gain, ||Cr||oo = K, find the optimal combination of the on/off rates 
and the length of the cascade that maximizes the signal amplitude, A, for any input R. 

To formulate this problem, first define the family of cascades that have the same internal 
gain K: 

C Kil := {C = (n,ai,.. .,a n ,p 1: .. .,/?„): = K}, 

Kpi ■ ■ ■ Pn 

For each input R, and each leak factor £, define the set of "optimal" cascades, that is, those 
cascades which exhibit maximal signal amplitude: 

C max (^,i?) := {C G C Kit : A(C,£,R) > A{C,£,R), for all C G C K , t }. 

Then define the function 



n 1 

<r (n,p\, . . . ,P n ) := ^2 -m 

i=l 



and observe that it satisfies 



a(C, £, R) = y 1 + <7 (n, &,...,&,) + q(R). 

Finally, define the set of cascades that minimize a® over the family C K ( : 

C*(£,R) := {C GC M : a (n,p u . ..,#»)< a (n',(3[, ...,/£), for all C G C K , t }. 

Our first result states that in fact the sets C*(£, R) and C max (£, R) are equal, or in other words, 
that an optimal cascade will simultaneously maximize the signal amplitude and minimize the 
signal duration. 

Lemma 4.1 In the notation defined above, C max (£, R) = C*(£,R), for all inputs R G U and 
leak factors £ > 0. 

Proof. Fix any £ > 0, and any R G U. Recall the notation C = (n, a±, . . . , a n , (3±, . . . , /?„). 
Given any C,C G C Kyt , the following equivalences hold: 

o- (n, (3i,..., (i n ) < cro(n', #,...,/%), 

^+a (n,(3 1 ,...,p n )+q{R) < ^ + a (n', /3[, . . . ,/3' n ) + q(R), 
a{C,£,R) <a(C',£,R) (17) 



and also 



a(C,£,R) <a(C',£,R) e> > fv||ff| ' 2 



a(C,£,R) ~ a(C',£,R) 
& A(C,£,R)>A(C',£,R). (18) 
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Therefore, (fT7|) and (JT8|) imply that, for any two cascades C,C G C Ki( , 

a (nJ 1 ,...^ n )<ao{n'j' 1 ,...,^ n ) O A(C,£,R) > A(C'J,R). (19) 

To show that C*(^, R) is contained in C max (£, R), pick any C G C*(£, ii). Then 

cr (n,/?i ; . . . ,/3„) < <T (n',/^, . . . ,/3' n ), for all C' G C K>1 . 

By (UU), this is equivalent to A{C,l,R) > A(C',£,R), for all C G C Kl<> and so C G C max (^,i?). 

Conversely, we need to show that C max (^, R) is contained in C*(£, R). So, pick any C G 
C max (£, ii). It satisfies: 

A(C,£,R) > A(C',£,R), for all C G C Ki< . 

Again by (|T9|) . this is equivalent to a"o(n, /3i, . . . , < ao(n' , (3[,..., (3' n ) for all C G C A -, £ . We 
conclude that C G C*(£, i?), as we wanted to show. I 

An immediate conclusion from Lemma 14. II is that, 

maximize A(C, £, R) over C Kt i 44> minimize o"o(n, f3\, . . . , j3 n ) over C Kit , 

so that, for any fixed internal gain, maximal amplitude is achieved simultaneously with minimal 
signal duration. This is consistent with the notion that the most efficient cascade would respond 
with sharp (high-peaked and fast) output signals. In the limit, this notion can be regarded as 
an "instantaneous response" (a « 0) coupled with "infinite signal amplitude" (.A oo), which 
is, of course, not biologically viable. A realistic solution to problem (P) does exist, and is stated 
in Theorem ^ 

Since the signal duration depends only on the cascade length and the "off" rates, Pi, (besides 
the input term), we expect the "on" rates, ctj, to play a small role in maximizing the efficiency 
of the output response. So, for addressing the problem (P), we will consider two different 
assumptions on the available knowledge on the a^: either (a) all the a, have an equal, fixed 
value, a; or (b) the product of the «j is known, at some fixed ap. We will also assume that 
the "leak" factor I is fixed, since this parameter was added artificially and may be adjusted 
independently. 

Before stating the main Theorem, we need to introduce some notation. Define the function 
/ : (1, oo) — > (0, oo) to be 

m = e [( 1 + I) ln ( 1 + I)-I 

Some properties of this function are stated in Appendix^] For any real number M > 1, define 

|_MJ = largest integer less than or equal to M, 
\M~\ = least integer greater than M, 

which are also known as, respectively, the "floor" and "ceiling" functions of M. Observing that 
any real number M > 1, can be written as the sum of its integral and fractional parts: 

M= [M\+5 M , 

where 5m G [0, 1), define the function ^ : (— oo, oo) — > N, which is plotted in Figured 

1, M < 1 
t>(M) = { [M\, M > 1, and 8 M < f([M]) 
\M], M>1, and.5 M >/(LMj). 
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This is a step function where the "jump" between steps always occurs in an interval between 
two integers, say k and k + 1, at a point that depends on the number k. In particular, since the 
function / is strictly increasing and takes values in the interval (2 In 2— 1, 1/2) (see Appendix^)), 
it follows that in some cases only the fractional part of the number M affects the location of 
the "jump" discontinuity: 

< S M < 2 In 2 - 1, *(M) = |MJ , 
^<6 M <1, *(M) = \M], 

while for the other cases, 2 In 2 — 1 < 5m < 0.5, the choice depends also on the integral part of 
M. 




-1 1 2 3 4 5 6 10 20 30 40 50 60 70 

M KL 



Figure 4: Left: the function 'I'(M). Right: the function ^>{2\n.K£). Note that, for a given gain 
K and leak factor £, the optimal length is given by the integer platform corresponding to the 
product KI. 



Theorem 1 Let K > and £ > be fixed real numbers. Let C K ^ e be the set of all cascades (jSJ) 
with internal gain K, as defined above. Then 

1. For each fixed n = N £ N, the elements C = (N, a%, . . . , «at, f3±, . . . , [3n) £ R) 
satisfy = /3, for alH = 1, . . . , N, where 

P ~ V Ki ) ' 

2(a). Any element C G R) of the form C = (n,a, . . . ,a, (3±, . . . , f3 n ) satisfies 

i 

n = (2 In «) and /3 i = /3 = a(— 

\K£ 

2(b). Any element C £ C*(£,R) of the form C = (n, a±, . . . , a n , (3%, . . . , f3 n ) G C*(£,R) 
with ax ■ ■ ■ a n = ctp satisfies 
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Before presenting the proof of the Theorem, some remarks on the interpretation of points 1 
and 2(a), 2(b). The first part of the result is consistent with the observation that the ordering 
of the amplification or dampening single steps within the cascade does not influence the final 
output signal (also observed in [TU]). 

The second part of the Theorem shows that indefinitely increasing the cascade's length will 
not increase amplification. In fact, there is an optimal length for the cascade that provides 
both maximum signal amplitude and duration. A similar observation was mentioned in |1U| . 
and our Lemma l4.1l and Theorem ^ characterize the conditions for achieving this optimization. 
For each gain K and leak factor £, this optimal length is easily read out from Figure 0] For 
instance, a cascade with a 6 to 9- fold gain (and i = 1), is seen to have an optimal length of 4 
steps. Figure El illustrates Theorem Q for an 8- fold cascade gain. The Figure shows the results 
of two simulations of system (J5J), both with input R(t) = 5te~ 2t , but different lengths of the 
cascade. The various curves represent R, the concentrations of each kinase Xj, i = 1, . . . , n, and 
the output X n+ \. It is clear that, for the non-optimal n = 7, the output's amplitude decreases 
and the signal duration increases. Note that the output curve X% is more spread out across 
time and its maximum value is lower, than for the optimal case. 



x, 




time, t time, t 



Figure 5: Signal transduction cascade with R(t) = 5te 2t , with K = 8, £ = 1, aij = 1.2. The 
horizontal lines represent A. Left (optimal case): n = 4, ft ~ 0.714, i = 1, . . . ,n, A = 0.409, 
and cr = 3.059. Right: n = 7, ft « 0.892, i = 1, . . . , n, A = 0.389, and a = 3.210. 

Theorem ^ can be proved by successively solving the two optimization problems: 

(PI) For each fixed n, minimize do, over all possible choices of 0i, . . . , j3 n € (0,oo), subject to 
||G||oo = -K". 

(P2) Minimize <ro, over all possible choices of n G N and (3\,...,j3 n G (0, oo), subject to 
l|G||oo = K. 

Recall that we are assuming that either (a) all the have an equal, fixed value, a; or (b) the 
product of the a% is known, at some fixed ap. The solution of (PI) is equal for both cases, but 
the solution of (P2) is slightly different for (a) or (b). Thus, problem (PI) is part 1 and (P2) is 
the part 2 of the Theorem. As we will see, the solution of (PI) greatly simplifies the proof of 
(P2). 
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4.1 Solving (PI): proof of part 1 of Theorem [T] 

Given a cascade of length n, this problem consists of finding a set of n parameters Pi, . . . ,j3 n 
such that the function uq attains a minimum value at Pi, i = 1, . . . , n, i.e., 

11 111 1 



ft Pi PU ~ ft ft ft-! 

for every Pi, . . . , p n such that ||G||oo = K- 

II^IU = 7"? °7T = K K£pi---0n-ai---a n = 0. 

I Px-'Pn 

For simplicity, rescale the values to -Bj = 1//3?, and observe that 

1 (P.^Pn? ^'"^ 



B 1 ---B n v ™' V « 
Then, the problem consists of minimizing the function: 

Q 



F(B l ,...,B n _ l ) = B 1 + --- + B n _ l 



Bi - - B n _x 

over all possible choices of Bi > 0, z = 1, . . . , n — 1, where 

Kl 



a\ ■ ■ ■ On 

In Appendix El we show that the solution to this optimization problem is 

Bi = Qn, i = 1, . . . ,n- 1, 

which also implies: 

Q i 

B n = „-l = Q n • 

So, the choice of the "off" rate constants that minimizes <jq is to have Pi = Pi = ■ ■ ■ = p n = P, 
with 

a I (oi\ ■ ■ ■ a n \ k 



as we wanted to show. 



4.2 Solving (P2): proof of part 2 of Theorem [T] 

To solve the more general problem, we first show how its statement can be simplified. Given 
the value of a (respectively, ap), suppose that we have found a solution of (P2), i.e., an integer 
n* and a set of constants /?*, i = 1, . . . , n* satisfying 

a (n*,P* 1 ,...,P^)<a (n,p 1 ,...,P n ) (20) 

for any other cascade C = (n, a\, . . . , a n ,Pi, . . . , p n ) with ctj = a, i = 1, . . . , n (respectively, 
«i---a n = ftp). 
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We have already showed that 

ao(n*J*,...J*)<a (n*,f3l...,P**) (21) 

with 

-„ ( a x ■ ■ ■ a n * - 



2 



Kt 

and we know this choice yields the unique minimum of <jq for a fixed length n. So, it follows 
that the solution of (P2) must also satisfy 

/3*=P*, i = l,...,n*. 

This observation allows us to simplify the statement of problem (P2), and look only for solutions 
where all /3j's are equal: 

(P2)' Minimize a (n,f3, .. .,/?) = n/(3 2 , over n G N and (3 G (0,oo), subject to {a/[3) n = Kt. 
From the constraint ||G||oo = K we have 

case 2(a): = K£ & (3 = a Q^V ^o(n, p(n)) = ^ n {Kl)$. 

case2(b): ^ = Kt o => *o(n,/3(n)) = n 

In either case, to solve the problem, it is enough to minimize the function ln[o"o(ra, /3(n))]: 

Fin, M) = mrt + — M 
n 

over n G N, where M is a positive constant with value either 

M = 2 In K£, for case 2(a) (22) 

M = 2 In — , for case 2(b). (23) 
ap 

For a fixed M, let the minimizer of F(n, M) over n G N be 

n*(M) := {n en : F(n, M) < F(n', M), for every ri G N}, 
which is given by Lemma IC.ll (Appendix EJ) : 

n*(M) = *(M). 

Thus, for part 2(a) of the Theorem we have n = n*(21n Kt) = ^(21n Kt), and for part 2(b) 
we have n = \I>(2 In Kl/ap). The value (3 is given according to part 1. I 



As shown in the example of Figure El evaluation of ao at n* (M) yields a value which is 
actually quite close to the "true" oo(M, (3{M)). 
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10 15 20 25 



30 



3.6 



Figure 6: The function a (n, 0(n)), for K£ = 8 and a = 1.2; and the points (M,a (M, 0(M))) 
(circle) and (n* , ao(n* , 0(n*))) (star). 



5 Cascades with positive feedback 

In this Section we investigate the behavior of cascades under positive feedback. Assume that 
the last kinase, X n , also contributes to the activation of the first kinase: then the differential 
equation for X\ includes one more term and becomes 



dX 1 



a x R(t) + eX n -PiXt. 



We will assume that e is small enough: 

01 ■ ■ ■ n > e«2 ■■■a n . 

This guarantees that the cascade is stable with respect to small perturbations (that is, all the 
eigenvalues of the system's matrix A have negative real parts, see Section Q). 

We can compute the transfer function for the system with feedback (s > 0), just as we did 
in Section^ for a given cascade C = (n, ct\, . . . , a n , . . . , (5 n ) and any input R, and leak factor 
I > 0. We obtain: 



goo 



1 



cti • • • a Tl 



S + £ (s + • • • (s + Pn) - £OL2 ■ • • On 



(24) 



The infinity norm is again obtained for the case u> = (s = jcj) (see below, at the end of this 
Section): 



IGI 



ax ■ ■ ■ Oi n 



£ 01 ■ ■ ■ /3 n - EOL2 ■ ■ ■ a n 

Computing the signaling time (r), and the signal duration (a) and amplitude (A), we have 



£ 01 • • • Pn ~ £«2 • • • a r 



ds 



(25) 



s=0 



1 (0i ■ ■ ■ 0n) 2 [E£=i j? + £"2 • • • a n Y. i+j ^7 
£ 2 {0i---0n-ea 2 ---a n ) 2 



+ q(R), 



(26) 
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AUCJ,R) = (r \ m -, p|| 2 , (27) 

a(C, £,R) fa ■ ■ ■ (3, n -ea 2 ---a n 

Comparison of these quantities for the models with and without feedback leads to the following 
conclusions: 

1. the system with feedback exhibits higher internal gain; 

2. the system with feedback exhibits larger signaling time and signal duration r fb > r and 
0fb > a. 

So, for an arbitrary cascade, the existence of a positive feedback leads to a less sharp output 
signal: the signal transduction down the cascade takes a longer time, and the output signal has 
greater duration. 

On the other hand, the existence of feedback may be used to great advantage in the design 
of an optimal cascade: positive feedback (at a constant rate e) allows the cascade to be of 
shorter length and still have the same maximal amplitude/minimal duration. The results in 
Theorem Q are valid just as before, with suitable adjustements to some of the constants. Thus, 
we now have 

HA,, jy- o a (ai+eK£)a 2 ---a n 
Halloo = & Pl"-Pn = 



K£ 

and now, similarly to the proof in Section f4. 11 we set 



Qt 



K£ ■ 1 



{a\ + eK£)a2 ■ ■ ■ a n t 
which leads to the optimal value for = /3a,, i = 1, . . . , n 

r, _ f (ai + eK£)a 2 - ■ ■ a n \" 

^-y Ki ) ■ 

To find the optimal length of the cascade with feedback, note that 

l / Kl 

<7 (n, (3{n)) = nM£ , with M ft = 2 In 



(«i + eK£)a 2 ■ ■ ■ ot n 

Since < M, then also n^(M fb ) < n*(M). Therefore, we conclude that, for the cascade with 
feedback, 

3. for each fixed n, the value of the off rates that maximixes A (minimizes (To) over C*(^, R) 
is larger, (3^ > (3; 

4. the length of the cascade that maximixes A (minimizes do) over C*(£,R) is smaller, 

< n*. 

These results agree with what would be expected from a signaling pathway: indeed, the 
existence of positive feedback enhances the activation at each step, so a larger amount of the 
phosphorylated kinase will be produced; to keep this amount at a "weak" level, the phosphatases 
should increase their activity. On the other hand, since the amount of phosphorylated kinases 
increased, a smaller number of steps is required to produce the same signal amplitude as in the 
cascade with no feedback. 
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To compute the infinity norm ||G||ooj we nrs t n °te that the denominator of G(jto), which 
we will denote by den(G(jo;)), satisfies (by the triangle inequality): 



|den(G(>>))l > \J U + ?\ [b w + Ai| ■ ■ ■ + /3 n \ -ea 2 --- a n \. 



Also 



\ju + P n \---\ju + n \ = ^(u 2 +[3 2 1 )---(u 2 + Pl)>p 1 ---(3 n , 

for every w£l, where the equality holds if and only if oj = 0. Thus 

\den(G(ju))\ >l[Pi---(3 n -ea 2 ---a n ] = den(G(0)) > 0, 

where the last inequality follows from the assumption {3±- ■ ■ f3 n > ea 2 • • • ot n . Therefore, if the 
expression |den(G(j(j))| is minimized at uj = 0, then the function |G(ja>)| is maximized at u = 0, 
as we wanted to show. 



6 Signal delay 

The frequency domain approach for linear systems also provides an answer to certain problems 
involving delays and stability within a signaling cascade. For instance, if there is delay in 
transmitting the signal at any step along the cascade, then both the amplitude and the signal 
duration are not affected. Suppose that, at each step, there is a delay 5{ in the transmission of 
the signal. The differential equation becomes 



d_ 

It 



( *i(t) \ 

x 2 (t) 
\x n+1 (t)J 



A 



( X 1 (t-S 1 ) \ 
X 2 (t - S 2 ) 

\X n+ i(t - S n+ i) I 



+ B R. 



The Laplace transform of Xi(t — Si) is, from the properties listed in the Appendix, 



so that 



The transfer function becomes: 

1 



X 



i+l 



e- s5i X 



ai+1 -e-^Xi. 



s + (3i 



+i 



G{s) 



a± - ■ ■ a r 



■ ■ ■ e 



—sS. 



n+l 



But, for an imaginary number juj, |e 



-sSi 



1, so the norm ||G||oo is unchanged; and since 



1 when evaluated at s = 0, the signal duration and amplitude are also unchanged. This 



is not surprising, because in a linear system, delay simply causes a temporal translation of the 
signal, by a fixed amount, without affecting amplitudes. 
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7 Stability of cascades 



A signaling pathway is considered stable (see jlUp if small and random perturbations (those 
that do not consist of biologically relevant inputs) are not amplified. So, in the presence of small 
perturbations, the amount of phosphorylated kinases should not be allowed to grow very large, 
and should return to the stable state, with X, ~ 0, for all i = 1, . . . , n. Thus, the behavior of 
a signaling pathway in the absence of a relevant input always satisfies expression ©, that is, 
X; <C X tot j for each % = 1, . . . , n, and hence its stability may be established by analysis of the 
model ©. 

In the absence of an input (R(t) = 0), the point (X\, X 2 , . . . , X„ +1 ) = (0, 0, . . . , 0) := is 
an equilibrium point of system © , and the stability of this equilibrium determines the stability 
of the pathway. The equilibrium point is stable if all the eigenvalues of the matrix A have 
negative real parts. This is indeed the case for the system described by equations (JSJ). We know 
that, after a perturbation, the system will always return to 0. Moreover, we can estimate that 
a small perturbation will also generate a small response, since: 

ll^pertlb ^ K 1 1 -Rpert 1 1 2 j 

where k is a constant, equal to ||Gr||oo- 

For signaling cascades which exhibit a lower degree of kinase specificity, the problem of 
stability of the cascade (see ^Uj) becomes significant. If a kinase X affects both the downstream 
kinases and some upstream kinase, then the eigenvalues of A change, and stability is not 
guaranteed. Allowing for kinase non-specificity, a resulting matrix A could be of the form: 



A F 



-01 


£12 





• 


£ln 




a 2 


-02 





• 








£ 




-03 ■ 


• 








£ 


£ 




• 








£ 


£ 


£ 




-fin 














■ 


1 


l) 



and the signaling pathway is stable if 

all the eigenvalues of A e have negative real parts. 
Some relevant easy-to-compute examples are: 

• suppose that each kinase % is only allowed to activate its downstream kinases (e\ n = 0) 
£\2 = and e ^ 0); in this case it is not surprising that stability is not affected at all, 
because this situation corresponds to a lower triangular matrix, again with eigenvalues 

-A; 

• suppose that there exists feedback from the last activated kinase to the first kinase (ei n = 
£o> £12 = and e = 0); in this case, if £o satisfies 

a 2 ---a n 



01---0 n > £o«2 



a r 



O £ < 



then the eigenvalues of the matrix A £ all have negative real parts and the new cascade 
is stable. To prove this, suppose that there exists an eigenvalue of A e with positive real 
part, that is, a complex number A such that 



A = X re + j\i 



At-r ^ 0, 



(28) 
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and 



det(^ £ - XI) = (-l) n 4(A + Pi) ■ • • (A + fin) - e a 2 • • • a n ] = 0. 



(29) 



Then (|2*8|) implies |A + Pi\ > for i = 1, . . . , n and so 

|(A + /3i)---(A + / 9 n )-e Q!2 > \{X + Pi) ■ ■ ■ (X + P n )\ -e a 2 ---a n 

> Pi ■ ■ ■ P2 - £o«2 ■ • ■ OLn > 0, 

which contradicts equation (|2T))) . 

8 Conclusions 

By modeling weakly activated signal transduction cascades as linear systems and applying 
techniques from control systems theory one can identify the cascade's input-to-output transfer 
function and internal gain. Based on these properties, the concepts of signal duration, signaling 
time and signal amplitude may be defined in an intuitive and general form, for any input signal. 

Our analysis shows that, for linear cascades, signal amplitude and duration are, respectively, 
maximized and minimized simultaneously. So, a cascade can respond with signals that are 
both fast and exhibit high amplification. To achieve the highest amplification and the shortest 
duration response, the cascade should have all off rates equal to some value p. 

We also show that, for each fixed internal gain, there are finite values for the length of 
the cascade and the off constants that simultaneously maximize (resp., minimize) the signal 
amplitude (resp., signal duration). To achieve these optimal conditions, the optimal length 
should be given by the well defined step function This function ^ depends only on, and 
increases logarithmically with, the internal gain of the system. The off constants should all 
have the same value p. This optimal value P depends on the internal gain and the length of 
the system. 

In addition, our analysis shows that a positive feedback term on the cascade enhances the 
optimal design, by allowing the same signal amplitude and duration to be achieved with a 
shorter length and higher off rates. 

Finally, other issues, such as delay at each phosphorylation step, and the stability of the 
signaling pathway when there is a high degree of non-specificity among the kinases, are also 
naturally examined within this framework. The stability of the zero steady-state of the cascade 
with respect to small perturbations is established by checking that the eigenvalues of the matrix 
A all have negative real parts. 

A Properties of function f(k) 

The function / : (1, 00) — > (0, 00) 



f(k) = k 2 




has the following properties: 



1. / is strictly increasing; 



2. /(I) = 2 In 2- 1 0.386 and lim^^ f(k) = 1/2. 
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To prove property 1, notice that another expression for / is f(k) = k[(k + 1) In (k + l)/k — 1], 
and compute the first and second derivatives: 



dk 2 



-n , fc + 1 2A; + 1 
(2fc + 1) In — — + 



k 



k(k+l) 



It is clear that the second derivative is always positive, and hence the first derivative is strictly 
increasing. Since df/dk(l) = 3 In 2 — 1 > 0, it follows that the first derivative is also always 
positive and therefore the function / is strictly increasing. 

To prove property 2, the value /(l) is straightfoward, and for the limit as k — > oo, it is 
easier to consider x = 1/k and compute: 



lim f(k) = lim f(l/x) = lim 

k^oo x^O x^O 



(1 + x) ln(l + x) - x 



0' 



This indeterminacy can be solved by twice applying L'HopitaPs rule: 
...first time: 



ln(l + x) + (l + x) T ^-l 

> -, as x — > 

2x 



...second time: 



1 



l 

^r- -> -, as x -> 0. 
2 2' 



0.5r 




0.35. 



5 10 15 20 25 30 35 
k 



Figure 7: The function f(k). 



B Minimization of ctq 

Let Q be a positive real number and n > 2 an integer. Consider the function F : (0, oo)™" 1 
(0, oo) given by 



F(B\, . . . , 5 n _i) — i?i + • • • + S n _i + 



Bi-B 



n-l 



Lemma B.l The choice of Bi > 0, i = 1, . . . , n— 1 that minimizes the function F is: i?j = Q 1 ^, 
i = 1, . . . , n — 1. 
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Proof. First, we claim that the search for a point (Si, . . . , S n _i) where F is minimized can be 
constrained to the compact set: 



n 



l^i Q 1/n , nQ 1 ^ 



n n-1 

(30) 



To justify the upper bound of the interval, observe that 

F(Q l ^ n ,...,Q 1 / n ) = nQ 1 / n (31) 

and that, for any j = 1, . . . , n — 1, 

Bj > n Q 1/n => F(Si, . . . , S n _i) > n Q 1/n . (32) 

So, it is enough to look for a minimum of F in the region Bi < n Q 1 /™, i = 1, . . . , n— 1 (because 
inside this region there is at least one point - equation (]31(1 - where F has a lower value than 
anywhere outside of this region). 

To justify the lower bound, suppose that Bj < Q 1 ^, for some j = 1, . . . , n — 1. Then 
using the already established upper bounds 

F(Si, . . . , S n _0 > Q > i nl/n jr nl/n]n 2 = n Q 1/n , (33) 

Si • • • Bj ■ ■ ■ S n _i -J^Ql/n [ n Ql/n]« 2 

and similarly we conclude that it is enough to look for a minimum of F in the region Sj > 
J= T Q 1 / n ,foTi = l,...,n-l. 

The function F is continuous, in fact differentiable, in the compact set (J30|). and so F 
has (absolute) maximum and minimum values in this set. The maximum and minimum may 
be attained either at a critical point of F, or at the boundary points of (|3()|h Equations (|32|) 
and ()33|) show that the minimum is not attained at any of the boundary points. So the minimum 
will be attained at an interior point of the set (jliOjl . which must also be a critical point of F. 
The critical points of F are given by: 

dF _ 1 Q . _ 

dSj Sj Si • • • S n _i ' 

or, equivalently, 

Q 



where S* satisfies 



d\ - ■ ■ rfn-i 



1-- ^r=0 O B* = Q l ' n . 

b* Br 1 



Thus, there exists a unique critical point of -F, (S*,...,S*), which indeed belongs to the 
compact set (|3()|) . By the discussion above, this point must be the minimizer of F, as we 
wanted to show. I 
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C Minimization of F(n, M) 

For a fixed M, let the minimizer of F(n, M) over n G N be 

n*(M) := {n € N : F(n, M) < F(n', M), for every n' G N}, 
Lemma C.l Let M be any fixed real number. Then n*(M) = ^f(M). 

Proof. Since n*(M) is the minimizer of F(n, M) over the (positive) natural numbers, we start 
by computing the derivative of F(n, M): 

^(n,M) = --^M = ^ [n-M\. 

dn n n z n z 

We consider two distinct cases: 
• Case M < 

^(n,M)>0, foralln>l, 
dn 

so F(-,M) is a strictly increasing function and thus its minimizer over N is the smallest 
natural number, i.e., n*(M) = 1. 



• Case M > 



(n, M) = 44> n = M, 



and the derivative is negative for n < M and positive for n > M: in other words, the 
function F has indeed a minimum at n = M. However, in general, M is not an integer, 
so it cannot be a solution to our minimization problem. We should choose 

r i, m < i 

n*(M) = i L M J> M > 1, and F(|~M],M) > F([M\,M) 
[ \M], M > 1, and F(\M],M) < F([M\,M). 

Note that we pick n* = 1 whenever M < 1, since a cascade of length zero is meaningless. 

To further analyze this condition, observe that we can write, for M > 1, 

M = k + 5, [M\ = k, \M] = k + 1 

where A; > 1 is the integral part of M and (5 G [0, 1) is the fractional part of M. Now, the point 
S for which n* "jumps" from \_M\ to \M~\ can be found by setting 

= F(\M],M) - F([M\,M) = F(k + l,k + 6) - F(k,k + S) 

= ]n{k + l) + -^—(k + S)-\nk-hk + S). 
k + 1 k 

Simplifying this equation we obtain: 
fc+1 k+S 



In 



fc fc(fc + l) 



5 = jfe(jfe + 1) In 

~k + l 



4^ 



b = k A 
5 = k 2 



k + l 
k ~ 
k + l 



k 



In 



1 + 



1 



In 1 + 



= /(*)■ 
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Analysis of this function (see Appendix ^J), shows that / is positive and strictly increasing, so 
we have 

F(\M],M) - F([M\,M) >0 f(\M])-5>0. 
Therefore, we should choose 

[ 1, M< 1 

n* (M) = < [M\, M > 1, and 5 < f([M\) 
{ \M\, M > 1, and 5 > f([M\). 

This proves the Lemma. I 

D Dictionary: Laplace transforms and transfer functions 

For further details about these topics see, for instance, [3] and 0, jS], [T8] . 
Laplace transforms 

For a function X : (— oo, oo) — ► M. n (with |X(i)| < ce kt , for all t, for some positive constants c, 
A;), the Laplace transform is another function X : 1Z — » C™ defined as 



/oo 
e" s *X(t)dt 
-oo 



where 1Z C C is the region of convergence of the integral. For example, if X(t) = e _3< , for t > 
and -X"(i) = otherwise, then X(s) = l/(s + 3), and 7?. = {s = s re + jSj m : s re > —3} (j is the 
imaginary number \J — 1). 
Some of its properties are: 

1. For any constant matrix j4 E R nxn 

AX(s) = AX(s); 

2. The Laplace transform of the derivative of X is 



/oo 
e~ st X(t)dt = X(0) + s l(s) ; 
-oo 



AX 
It 



3. If X(t + (5) =: W(t) is a translation of X, then 

iy(s) = e" s5 X(s) ; 

4. The inverse Laplace transform is 

1 PSre+]00 

X{t) = / e st X(s)ds 

271"J J S re — 70O 



with s = s re + JSjm, where s re is chosen so that s re + jSi m is in the region of convergence 
1Z. 
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Transfer function 

Let A G M. nxn , B G R nxm and C G M pxn be matrices, and let X G R n , Y G M m , i? G M p , and 
consider the n-dimensional linear system with m inputs and p outputs: 

dX , s 

— = ^X + Sfl, (34) 

F = CI. (35) 

Applying the Laplace transform operator on both sides of the linear system (|34j) - (|35j) yields an 
algebraic equation relating the new functions X(s), Y(s) and R(s): 

sX(s) = AX{s) + B R(s) 
Y{s) = C X(s). 

Moreover, for every s for which the matrix si — A is invertible (/ is the identity matrix), 

(si — A) X(s) = B R(s) => X(s) = (si - A)- 1 B R(s) 
and thus, one can solve immediately for the output 

Y(s) = C(sI-A)- l B R(s). (36) 
The transfer function of the system (|34j) is 

G(s) := C(sI-A)- 1 B, 
and depends only on the internal structure of the system (i.e., A, B and C). 

Impulse response 

A useful case is that of the impulse response: 

R(t) = 5(t), => R(s) = 1 

and therefore: 

Y(s) = G(s) & Y(t) = G(t), 

so that the transfer function of the system is the output corresponding to a single pulse of 
input. 

The gain H^Hoo 
We have 

-i roc -i roo 

\\ni = ir \G(ju)R(ju)\ 2 cLj<—sup\G(3uj)\ 2 / \R(jco)\ 2 du 
which is equivalent to 

11*112 < ||G||oo||«||2 \\Y\\ 2 < ||G||oo||fl||2. 
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So, the infinity norm of the transfer function is an upper bound on the strength of the output. 

To see that it is indeed the least upper bound, see for instance |HJ: we can always choose a 
frequency ujq so that 

II^Uoo = \G(ju )\. 
In our case, this is ujq = 0. Then choose a control such that 

= { 0, ftheiwise, (37) 

where e > and r should be such that R has unit 2-norm, for instance r = \J~kJe. For very 
small e > 0, |i?(ju;)| is zero, except on a very small neighborhood of loq = and we may 
approximate: 

1 f°° - 1 f £ 

— \G(ju)\ 2 \R(ju;)\ 2 dco « — / r 2 \G(j^)\ 2 dco 



e 



^\G(ju )\ 2 j r ,U 



—£ 



\G\ 



oo 



where the last equality follows from the definitions of loq and r. Therefore 



\Y\ 



1 f°° 

- I \G(ju)\ 2 \R(ju)\ 2 dco 

— oo 



2tt 



IGII 



As an example of an input that (approximately) satisfies Q37|) . consider R(t) = 2^sinei (for 
t > 0), the input plotted in Figure |3 Computation of the Laplace transform yields R(s) = 
^ [-7T — Arctan|] , where the function Arctan is the principal branch of the complex inverse 
tangent function. It can be shown that, for sufficiently small e, the function R approximately 
satisfies condition (|37|). except at the discontinuity points u> = ±e. 

More generally, in a system with m inputs and p outputs, one defines the internal gain of 
the system ||G||oo as the induced C 2 operator norm of the map from the inputs to the outputs. 
It is possible to prove that 

\\G\\ O0 = sap5[G(ju)] 
where 9 denotes the largest singular value of the matrix G(jco). 
Stability of the transfer function 

As remarked above, expression ()36|) is valid if and only if the matrix si — A is invertible, or 
equivalently 

s ^ A, for every eigenvalue, A, of A. 

If A m = max{Re(A) : A is an eigenvalue of A}, then the region of definition of the transfer 
function is included in the set 1Z = {s = s re + jSi m : s re > A m }. 

If all the eigenvalues of the matrix A have negative real parts, then the transfer function is 
said to be stable. This is case for the matrix of the signaling cascade whose eigenvalues are: 
— j3i, . . . , — /3 n , so the transfer function G(s) is stable and well defined on 1Z = {s = s re + jSi m : 
s re > - rain Pi}. 
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